Scalable Superconducting Architecture for Adiabatic Quantum Computation 
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A scalable superconducting architecture for adiabatic quantum computers is proposed. The ar- 
chitecture is based on time-independent, nearest-neighbor interqubit couplings: it can handle any 
problem in the class NP even in the presence of measurement errors, noise, and decoherence. The 
implementation of this architecture with superconducting persistent-current qubits and the natural 
robustness of such an implementation to manufacturing imprecision and decoherence are discussed. 

PACS numbers: 03.67.Lx, 85.25.Cp 



Adiabatic quantum computation is an approach to 
solving computational problems of the complexity class 
NP |3| via energy minimization. In particular, by ex- 
ploiting the ability of coherent quantum systems to fol- 
low adiabatically the ground state of a slowly changing 
Hamiltonian, it aims to bypass the many separated lo- 
cal minima that occur in difficult minimization problems. 
Adiabatic quantum computation is of theoretical interest 
because it provides a straightforward, non-oracular way 
to pose class NP problems on a quantum computer. To 
date, most research on it focuses on ascertaining its time 
complexity |3j S S IE 0- However, it is also of prac- 
tical interest because encoding a quantum computation 
in a single eigenstate, the ground state, offers intrinsic 
protection against dephasing and dissipation 0, 01 . 

In this Letter, we present a scalable superconducting 
architecture for adiabatic quantum computation that can 
handle any class NP problem. It requires neither effi- 
cient qubit measurements, nor interqubit couplings be- 
yond nearest neighbors, nor couplings that vary during 
the course of the computation. We also discuss how to 
implement this architecture specifically with supercon- 
ducting persistent-current (PC) qubits which con- 
stitute a promising approach to lithographable solid-state 
qubits. We show that the proposed architecture is robust 
against manufacturing imprecision, and we estimate the 
maximum size problem the architecture could support 
if it were implemented with existing PC qubits at dilu- 
tion refrigerator temperatures of 10 mK. This maximum 
stems from the condition that the environment must not 
excite the computer from its ground state. A simple 
Boltzmann factor argument implies that a 10 mK tem- 
perature limits the PC qubit architecture to NP problem 
instances of 0(50) bits. However, advances in cryogen- 
ics or in fabricating PC qubits from higher-Tc materials 
could raise this limit by 1 or 2 orders of magnitude. 



Background: In adiabatic quantum computation, one 
encodes the answer(s) to a hard constrained minimization 
problem in the ground state(s) of a suitable Hamiltonian 
H whose local couplings ensure its ground state(s) sat- 
isfy the problem's constraints. One initiates the system 
in the ground state of another Hamiltonian i?o chosen so 
its ground state is quickly reachable simply by cooling. 
One then adiabatically deforms Hq into H by applying a 
time-dependent Hamiltonian H{t) such that H{Q) — Hq 
and H{T) = H. The adiabatic approximation holds as 
long H{t) possesses at all times a spectral gap A{t) be- 
tween its instantaneous ground \'ipo{t)) and excited states 



\ipn{t)) such that 
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of whether adiabatic quantum computation has polyno- 
mial or exponential time complexity is thus determined 
by whether the minimum gap Aniin(?^) shrinks polyno- 
mially or exponentially in the number of qubits n. 

It is still unknown what speedup adiabatic quantum 
computation offers in general over classical energy min- 
imization algorithms. Numerical investigations of hard 
instances of NP-complete problems 0, suggest that 
Amin('T') = 0{n~^), at least typically, and thus adia- 
batic quantum computation may provide for all prac- 
tical purposes an exponential speedup on NP-complete 
problems over all known classical algorithms. However, 
quantitatively bounding the minimum gaps for adiabatic 
algorithms for NP-complete problem instances is appar- 
ently as hard as solving the actual NP-complete instances 
themselves. Thus, these numerical investigations have 
been confined to problems with < 30 qubits. Ref. 
constructs a problem that takes exponential time for 
a sirnple version of the adiabatic algorithm. However, 
Ref. shows that alternative versions of the adiabatic 
algorithm solve the problem of ^ in polynomial time, 
and presents a problem for which simulated annealing 
provably takes exponential time yet adiabatic quantum 



computation takes only polynomial time. Definitively 
establishing the efficiency of the adiabatic algorithm for 
large problem instances may thus have to await the ac- 
tual construction of adiabatic quantum computers of the 
sort described here. 

Layout Requirements: By definition, a scalable pro- 
grammable architecture for adiabatic quantum comput- 
ing possesses some regular geometrical layout of switch- 
able pairwise couplings on m nodes that can efficiently 
encode problem Hamiltonians corresponding to any in- 
stance of some desired class of energy minimization prob- 
lems up to n < m bits in size. To have the maximum 
flexibility in the problems one can pose and to exploit the 
maximum potential power of adiabatic quantum compu- 
tation, the binary constraint problem the architecture 
naturally poses should be NP-complete, and the archi- 
tecture should readily pose instances of it that are truly 
difhcult for all known classical heuristics. Additionally, 
present experimental constraints make it highly desir- 
able to have couplings that neither extend beyond near- 
est neighbors nor vary in time beyond possibly switching 
between a full strength "on" state and a much reduced 
strength "off" state before computation begins so as to 
allow programming of the desired problem. 

An architecture mee ting these requirements naturally 
follows from the fact that calculating the ground 
state of an antiferromagnetically coupled Ising model in 
a uniform magnetic field is isomorphic to solving the 
NP-complete graph theory problem Max Independent 
Set (MIS), which is the problem of finding for a graph 
G ~ {V, E) the largest subset S of the vertices V such 
that no two members of S are joined by an edge from E. 
Interqubit couplings may be laid on a single plane and 
kept to a modest number since MIS remains NP-complete 
even in the ostensibly simple, topologically uniform case 
of degree-3 planar graphs, i.e., graphs that can be drawn 
in a plane without any edges crossing and in which every 
vertex is connected to exactly 3 others [l^ . 

The isomorphism takes a particularly simple form in 
this case: the MISs of a degree-3 planar graph G — (V, E) 
are the ground states of an Ising model with spins on ver- 
tices V , equal strength antiferromagnetic couplings along 
some desired axis h on edges E, placed in a uniform ap- 
plied magnetic field also along fi: 

In regard to the requirement of posing truly hard prob- 
lem instances, note that the most efficient classical ap- 
proximation algorithm known for MIS restricted to pla- 
nar graphs has a cost that grows exponentially in the 
desired accuracy. Specifically, the cost to obtain an ap- 
proximate MIS of a planar graph with n vertices that 
has a cardinality which is at least of the true MIS's 
cardinality is 0{8^kn), and is thus impractical once one 
desires > 90% accuracy (k > 9) [l3| . 
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FIG. 1: Schematic of tlie embedding of the simplest degree-3 
planar graph (left) into the architecture (right). The architec- 
ture is built on a triangular lattice with qubits at its vertices 
(schematically denoted by filled circles if possessing a Zeeman 
term and thus denoting a graph vertex, open circles if not and 
thus denoting a dummy vertex that serves only to propagate 
a couphng) and switchable nearest-neighbor couplings on its 
edges (schematically denoted by black double lines if antifer- 
romagnetic and switched on, black single lines if ferromagnetic 
and switched on, and gray broken lines if switched off). This 
allows an arbitrary degree-3 planar graph to be embedded ef- 
ficiently within the lattice, as can be seen in how the Qubit 
#3 of the graph is ferromagnetically coupled via dummies to 
a partner #3' so that is a nearest neighbor of #2. Inset: 
To compensate for measurement errors, the logical qubits in 
the architecture (denoted by the larger circles) are actually 
ferromagnetically coupled clusters of physical qubits (smaller 
circles). For example, as depicted above, one could ferromag- 
netically couple 6 dummies to each physical computational 
qubit and thus essentially compensate for measurement er- 
rors by a classical 7-bit repetition code. The ferromagnetic 
couplings in such redundant clusters need not be switchable. 

The requirements of nearest neighbor couplings, max- 
imum uniformity, and time-independent control can be 
met via the layout depicted in Fig. 1: a triangular lattice 
with qubits at its vertices and nearest-neighbor ferro- or 
antiferromagnetic couplings on its edges. To allow pro- 
grammability, that is, to embed an arbitrary degree-3 
planar graph in this triangular lattice, each coupling is 
switchable and moreover the magnetic field on each qubit 
is individually controllable so as to allow one to create 
"dummy" qubits that possess no single qubit Hamilto- 
nian and act simply to propagate a coupling. (Triangular 
lattices make it especially easy to embed degree-3 planar 
graphs, but lattices with a lower node degree, such as 
square or hexagonal lattices, would also suffice.) 

Implementation with PC Qubits: Fig. 2A depicts the 
specific PC qubit circuit we shall consider here. Ref. 
explains the rationale underlying its design and its canon- 
ical quantization. For our purposes, the pertinent result 
of Ref. 10] is that in the regime of sufficiently low temper- 
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FIG. 2: (A) Circuit schematic for PC qubit, after Ref. [T3|- 
(B) Schematic (not to scale) of two qubits coupled by induc- 
tive coils connected by a bus possessing a switch such as a 
SQUID (depicted) 0] or a JoFET O. Qubits are separated 
so that free space inductive coupling is negligible, and cou- 
pling via the inductive coils dominates. If the coils trap zero 
net flux, if Mtq denotes the mutual inductance between the 
flux transformer and a loop of the qubit, and if Lt denotes 
the self inductance of the flux transformer, then the effective 
mutual inductance between the qubits Mqq ~ M^q /Lt and 
the effective qubit self inductance Lq is largely suppressed. 



atures (T < lO^^Ej « 30 mK for our parameters) and 
frustrations that make all the minima of the junctions' 
potential energies nearly equal (0.485$o < f^°*- + lf^°P < 
0.515$o); the circuit is accurately described by a trun- 
cated 2-lcvel qubit Hamiltonian that has a simple, readily 
tunable form. Specifically, numerical analysis shows for 
realistic parameters of Ej/Ec — 80,/3 — 0.8,7 — 0.02 
and frustrated at /t°p(bot) ^ [0.330 + ^*°p(''°*)] $o yields 
an effective qubit Hamiltonian Hq: 

(2) 



where 



K, (5*°P,(5^°') 



-0.025 + 3.8^^°* + 2.05*°P (3) 
0.0049 - 1.2(5''°* - 0.81(5*°P (4) 



and the "z-basis" is the basis of classical counter-rotating 
current states, {|0), |0)}, while the "x-basis" is the basis 
of real symmetric and antisymmetric linear combinations 
of the classical states, {|0)+ |0), |0)- |0)}. 

Beyond this regime where Equ. (5) is valid (i.e., 
^top _|_ ijbotj > 0.015), the circuit overwhelmingly favors 
a single circulation direction for its current, meaning in 
the reduced 2-level picture that the component dom- 
inates over the ax by well over a factor of 20. Such a 
CTj component can dominate over all couplings, naturally 
providing a useable starting Hamiltonian with a ground 



state that is reachable by simple cooling despite the pres- 
ence of couplings that are always switched on. 

Qubits are coupled inductively as in Fig. 2B. Switch- 
ing of couplings between a full strength "on" and a much 
reduced strength "off" could be accomplished, for ex- 
ample, magnetically via DC SQUIDs Q or electrostat- 
ically via JoFETs (Josephson Field-Effect Transistors) 
[l^ . Now consider an inductively coupled pair of iden- 
tical PC qubits with identical applied frustration offsets 
= (5'°P,(5!;'(°2*) = 5^°^- Moreover, for simplicity, let 
the mutual inductances from each qubit to either loop of 
the other be identical, i*2^2i)'*^ = Inductive effects 
are calculated perturbatively to lowest order by equating 
the qubit circulating current to the mean current trav- 
elling through the DC SQUID in the circuit, which the 
previously cited single qubit numerical simulations show 
is roughly ±1.4/^ cos[7r(0.330 + 5% 

Within this approximation, the basic building block of 
our desired problem Hamiltonian of a Max-Independent- 
Set-encoding antiferromagnetic Ising model in a uni- 
form field, Equ. (3), is achieved at an operating point 
((5^ = -0.0124,(5^ = 0.0200) with Ej/{MI^) = 90: 



He 
Ej 



0.029 (tTcl + (Jc2 + CFcIO^ct) 



(5) 



where c is the vector rotated 19° clockwise from the z- 
axis in the xz-plane of the Bloch sphere. 

Realistic parameters for a Nb PC qubit are Ej = 
0.60 THz and Ic = 1.2/iA. Therefore, the desired ratio 
Ej/{MI^) — 90 requires a mutual inductance of M = 3.1 
pH, which is also realistic. 

We now turn to the problem of constructing dummy 
qubits that only propagate ferromagnetic couplings and 
do not possess single qubit terms. Ideally, these ferro- 
magnetic couplings would be proportional to —crci<7c2 
and thus commute with Hq- However, if we constrain 
our inductive couplings to have only 2 settings: "on" 
where M = 3.1 pH and "off" where A/ « pH, then the 
operating point ((5^ = -0.0171, = 0.0152) that causes 
the single qubit terms to vanish will not yield a coupling 
Hd oc —ac\(Jt2- However, this poses no practical problem 
as the actual coupling is very close to —ac\Oc2'- 

^ = -(3.6 X 10-3) a^^a^^ + (5.4 x 10-5)a,ia,2 (6) 

where d is the vector rotated 16° clockwise {i.e., 3° less 
than c) from the z-axis in the xz-plane of the Bloch 
sphere. (It probably unwise at present to expend any ef- 
fort toward making Hjj closer to —(Tci<Tc2 given the likely 
level of accuracy of our formulas for the Hamiltonian.) 

Similarly, the Hamiltonian Hcd for a computational 
qubit inductively coupled to a dummy qubit is not in the 
ideal form — o'ciO'j27 but again it is tolerably close: 



H, 



CD 



Ej 



0.029aei - 0.013crgicre2 



(7) 
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where e is the vector rotated 24° clockwise {i.e., 8° more 
than d) from the z-axis in the Bloch sphere's xz-plane. 

The adiabatic computation is performed simply by 
slowly bringing (5*°^, J*^"') on all the computational 
qubits from any convenient point with i5*°p + ^5^°^ > 
0.015 to the operating point (5^^,(5^°*) while keeping all 
dummy qubits fixed at their operating point {6*'^^ , ^^^)- 

The intrinsic robustness of adiabatic quantum compu- 
tation versus environmental noise has been demonstrated 
both numerically and theoretically 0. However, mea- 
surement error is presently a critical concern with PC 
qubits because their most conveniently measured char- 
acteristic is the flux created by their persistent currents. 
This limits them to be measured in what we have called 
the 2;-basis despite the fact that it often will be necessary 
to work it a different basis for computation. Specifically 
in the case of this architecture, the computational ba- 
sis is rotated 19° away from the z-axis. Therefore, any 
measurement scheme for this architecture based on mea- 
suring a qubit's magnetic flux will have an error proba- 
bility of at least sin^(19°) = 0.11. Moreover, a special 
concern arises in adiabatic computation based on frus- 
trated Ising models since such systems generically have 
highly degenerate ground states. Simple repetition of the 
algorithm will thus generically yield a different ground 
state with each measurement. Measurement errors can- 
not be corrected by averaging such uncorrelated data to- 
gether. It is therefore necessary to program correlated re- 
dundancy into the architecture by having dummy qubits 
ferromagnetically coupled to each computational qubit. 
Thus, when measurement collapses the computer's state, 
one will obtain multiple copies of one valid solution, and 
then averaging the data from measuring all these added 
dummy qubits will compensate for measurement errors 
via a classical repetition code. The couplings between 
these redundant dummies and the computational qubits 
need not be switchable, and can be designed in such a 
way that the addition of redundancy does not decrease 
the minimum gap Amin. 

The key constraint on the number of qubits in an adi- 
abatic quantum computer is that the environment must 
not excite the computer out of its ground state. Conser- 
vatively, this imposes a limit rimax on the number of logi- 
cal qubits such that Aniin(nmax) > kT, the environment's 
average thermal energy. As the form of A^in{n) is still 
an open problem, no exact answer can be given presently. 
However, as cited previously, there are indications that 
A^in{n) = O(n-i), at least typically [1I1I3. If this is 
true asymptotically, then the maximum number of qubits 
given an operating temperature T is 0[Ai/{kT)] where 
Ai is the energy gap of the problem Hamiltonian for a 
single computational qubit. For the previously cited pa- 
rameters for existing Nb PC qubits, Ai = 18 GHz. As di- 
lution refrigerators can bring the electron temperature to 
15-20 mK, the above assumptions imply rimax = 0(50). 
Such a value essentially is the maximum possible at dilu- 



tion refrigerator temperatures with Nb PC qubits of any 
possible design for it saturates the limit on Ai set by the 
Nb's critical temperature Ai = O(O.lfcTc) = 0(0.9 K). 
Advances such as building Josephson junctions out of 
higher-Tc materials and/or achieving electronic tempera- 
tures in superconductors of 0(1 mK) could increase this 
conservative estimate of n^ax to hundreds or thousands. 

Finally, the architecture is robust versus manufactur- 
ing imprecision for three reasons. First, any undesired 
term in the Hamiltonian that is a product of Pauli oper- 
ators will couple a problem's solution to only one other 
excited state. Perturbation theory therefore implies that 
for given energy scalar k [see Equ. (5)], tolerance for in- 
accuracies in the couplings decreases only linearly in n. 
Second, since all PC qubits have two controllable pa- 
rameters, (5*°P and 6^°^, certain fabrication errors can be 
compensated by individually calibrating ^*°p and 5^°* for 
each qubit. Third, if some qubits completely malfunc- 
tion, the redundancy in the triangular grid layout allows 
one to steer coupling chains around them. 
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